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SYNOPSIS 


P^cision Making is the process of selecting an 
alternative among the available ones. When there arc 
more tha.n one criteria for the selection, the problom is 
called multi -criteria problem. Most real life problems 
fall under this category. There are various methods to 
solve this class of problems. All those methods can be 
classified into four m.ajor catagories depending on the 
preference information provided by the decision maker end 
the stage at which it is provided. 

In this thesis wo develop an algorithm for bicriterion 

posterior 

linear programs. The algorithm requires/ articulation of 
information. The algorithm generates all officient extreme 
points in the objective space. Almost all the existing 
algorithms for these class of problems generate efficient 
points in the d-^cision spnee. However, the a-lgorithm proposed 
by Aneja and Hair [1] concentrates on the generation of 
efficient points in the objective spoce. Hence the proposed 
algorithm is compered with algorithm proposed by Aneja and 
Hair. For both the algorithms computer codes in FORTRAH 
l!\r were developed and implemented on DEC 1090, Computer 
system. 
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j?or coTnputxtion'-'.l simplicity, Ixlanfeci transportation 
problems of various sizes ’■"Gro developed. These v^ero solved 
as linear programs. The size of the line-ir program varied 
from problems with 7 constraints 16 variable to la constr'^'ints 
49 variables. 

The execution times for the gener-tion of one extreme 
efficient point for both the algorithms were compared for 
2i randomly generated problems of 4 different sizes. It was 
observed that the proposed algorithms is computational! 3 ^ highly 
efficient as compared to Anoja and Nair's algorithm. For the 
smallest size problem, i.e., 4x4 ba,lanced transportation 
problem, the a,vorage execution time (average b.ased on 
6 problems) per extreme efficient point w as I /6 that of the 
Aneja and Fair's algorithm. For biggoi sizod problems, 
the computo.t ional efficiency was even still better. 


/ 



CHAPTBE I 


IITTROJlUCTIOIsr 

Decision .naking la 3 , process of s^'locting a possible 
course of action among thv^ various ava,il 0 .ble ones. There 
may bj one or more crit .ria uhich xjill h-'lp justify the 
soloction. If thcT'- arc more than one criterion then the 
process is called Multipl'- Crit'.-ria Decision Making (MCDM). 
Formally the cri+cris art- non-commensurable and conflicting. 
Of la 1.0 ther-o has boon considerable research activity in 
th-. area of multi-critrT ia d-cision making sine-,.- most of the 
real life problems fall under this category. In the MCDM 
problems, tho best nlternativo is to be selected keeping 
in thr- O') ncidor.a tion, the various con,str=’ints .and criteri.a. 

In the process of choosing the best alternative, the d 'cision 
maker is raouirod to give some preference information. The 
stage at which it is given is the basis on which tho MCDM 
problems can bo classified. Th classific ition of this 
problems is given in Tabic 1.1. 

'i'he methods for whicn no articulation of pri^f er>.-nce 
information is given, assura" that the decision maker will 
be abl.”' to accept the given solution. This asa'omption is 
not well justified. The methods whore the decision maker 
is required to give his proferonco information apriori, 
suffer from the disadvantage that the decision maker gives 
this in an inf or’T''ation vacuum. The methods which are based 
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on ordinal prcft-rencc information given apriori depend 
heavily on this and the moment the decision maker finds that 
tho information given at an earlier time does not hold a,ny 
morc'j, tho problem is to bo solved afresh using the new 
information. 

In the third method, which is classified as interactive 
mf-^thod, the decision mak-^r is part of the solution process, 
Thor, is no nc.,d of a priori preforonce inforTr -ion arti- 
cula.tion. The majoi' advantage of this is that decision 
maker learns about the problem as ho proceeds. Since decision 
maker is part of tho solution process, the prospects of the 
obtained solution being implemented, is bright. But thw 
solution depends upon the accur'^cy of the- local preforonce 
tho decision maker can indicate. This is one of hie disadvan- 
tages, these methods suffer from. Another disadvantage is 
that much more effort is required on the part of tho decision 
maker than other methods. 

The last category of methods, listed in Table 1.1 are 
the methods for which a posteriori artifulation of prefe- 
rence information is given. Those methods hi^vo an advo.ntag-. 
over other methods, because no assumption regarding tho 
decision maker’s utility function is required. 

The present thesis deals with this class of multi- 
criteria problems. Solving thu mul.ti-crit-ria problem 
requires the generation of al.l non—dorninated (efficient) 

points. 
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Tho first generaliz-^tion of tho multi-critc.ria 
problnm is the hicriterion problems. If all the cons- 
traints, and both the criteria arc linear, than thr 
problem is called bicriteria linear program. The 
importance of the bicritcrion problems can hnot be over 
empha.sizod. Many situations co\iLd be formulated as bi cri- 
teria linear progra,ms. For instance in a production 
system one would like to minimize the cost of pi'-oduction 

and maximize the pin duct ion. 

A survey of the litora.ture indic'itos that most of 
the investigators [2, 4 » 9,10] ha.ve approached the bicriteria 
problem through developing algorithms for generating all 
th" efficient points in the d-cision space. These algori- 
thms involve the consideration of n dimensional space, 
where n is the number of decision variables. However, 
the dimensionality of the solution space can be considerably 
reduced if instead of the decision spac€.‘, the objective 
space is considered. In the present thesis, an algorithm 
whi.ch generates the extreme efficient points in the objective 
space is developed. The perfoimance of this algorithm 
is compared with Anoja and IT.iir's algorithm [l] the 
only other algorithm reported in the literature which 
exploits objective space rather than oocision space for 
the generation of extreme efficient points. A s-t of 24 
randomly generated balanced transportation problems of 
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dirff;:'cnt sizes are considrrud for comparing the two algo- 
rithms. For "both tho algorithms ^ those transportation 
prohloms are structured as linear programmes. 

The thesis is organized in four chapters. Chapter II 
deals with problon statement, its formulation and the 
relo'vant literature survey. The proposed algorithm is 
presented in Chapter III. In Chapter lY, a comparison 
of the proposed algorithm with Anoja and Nair's algorithm 
is presented. Listing of o.3.1 the computoj* programmes 
developed for this work ?r'P givon in Appendices A, B and C. 



CHAPTER II 


PPOBLEM STATEiyiENT AND LITBPA.TUEB SURVEY 


2.1 Problem Statement ; 

Bicriterion linear Programs (BLP) may be stated 
as follows ; 

Consider a Vector Minimum Problem (VTCP) 

^Cx, Dx 
s.t. Ax 

X 7/0 

where x, C and D are n-dimensional vectors, A is a (m x n) 
matrix and b is m-dimensional vector. C and D are coefficient 
of the two criterion functions. A represents the technolo- 
gical coefficients of the activities and b is the require- 
ment vector. AX ^ b represents the constraints of the 
linear programme (IP). This statement is made without loss 
of generality, as maximisation problems and greater than 
or equal to constraints can be readily incorporated using 
standard techniques of IP. 

let us denote the set of feasible solutions by 
S and its mapping by the two objectives into objective or 
criterion space by T. Y is called the pay off set and 

is defined as 7 

Y = zj) I Si = Oi, Z 2 = Dx for some x e 
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A solution is called superior if it minimizes 
both the objectives simultaneously. In most real life 
situations superior solutions do not exist. The objectives 
are non-commensurable. Under those assumptions the 
decision maker is interested in efficient solutions or 
non-dominated solutions. For clarity we define efficient 
solution before proceeding further. 

Definition ; 

A feasible solution x° is said to be efficient 
(non— dominated ) if Gx Cx'^ for some other feasible x 
implies that Dx°< Dx and vice versa. 

Essentially it means that no other feasible solu- 
tion is better with respect to both the criterion. Also 
by implication x e S is dominated if there exists atleast 
one feasible solution x° such that Cx° Cx and Dx°< Dx. 
Thus, it is always wise to restrict ourselves to non- 
dominated solution. The solution to (l) can be interpreted 
as determination of the s .t of entire efficient points. 

2.2 Survey of Literaturej, 

Most of the research in this area lias been concen- 
trated on getting all the efficient points in the decision 
space. One of the early papers in this area is due to 
Geoff rion [4]. He proposed a scalar maximization approach 
for obtaining all the efficient solutions. He converted 
the original problem into the following proolem, 

I 
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Problem Pq.? 

Min. aCz + (l-a) Dx 

S . t . Ax 4 b 

X ^0 (2) 

where a <£ (0, 1). The problem P^ gives all efficient 
solutions if a is varied in the given range. He proves 
the necessary and sufficient conditions in the form of 
a theorem. Further, G-eoffrion suggests that his 

approach can be extended to multi objective programs with 
more than two criteria. Bacopoulas and Singer [2] adopted 
a constraint criteria approach to get all the efficient 
solution in the decision space. One criteria is converted 
into a constraint at certain level. By parametrically 
varying this level he proves that all efficient solution 
can be obtained. 

Steuer [9] has given an algorithm which is a modi- 
fied version of the simplex algorithm. In the column selec- 
tion rule stage, the efficiency of the next solution is 
tested by introducing each possible candidate colimns 
one by one. One LP is solved to test the efficiency or 
non domination of each point, They call this procedure as 
adjacent efficient basis algorithm, 

Yu and Zeleny [10] have given another modification 
to the simplex procedure. They r ef er their modified simplex 
procedure as adjacent non-dominated -basis approach. A 
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modified version of the column selection rule is used to 
examine non-domination, ilssentially they use the c . of 

J 

each objective to test non-dominance. 

Both Steuer and Zeleny and Yu concentrate on 
generation of efficient points in the decision space, 
which requires more computing and book keeping. Purther, 
they have addressed to multi criteria problem in general. 

Sadagopan [8] has suggested an interactive solution 
procedure for the bicritcria mathematica.! programs. 

Recently, Aneja and Hair [1] have developed an algorithm 
basically for bicriterion transportation problems. They 
claim that their algorithm is, applicable to bicriterion 
linear programs also. 

Their algorithm generates efficient extreme points 
in the objective space. The algorithm requires solving of 
(2K-5) linear programs if there are K ( '^ 2 ) efficient 
extreme points in the objective space. As the present work 
considers the same problem and provides algorithm as com- 
pared to Aneja and Hair’s approach, the details of their 
algorithm are presented in Chapter III. 

The proposed solution methodology involves the 
generation of efficient extreme points in the criteria space, 
instead of decision space. This was done for the following 
reasons. Firstly, the dimensions of the solution space get 
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reduced if the objective space is considered. Por instance 
in the case of bicriterion problems the objective space is 
two dimensional while the decision space is n-dimensional 
where n is the number of decision variables. Secondly, 
all extreme points in the criteria space correspond to 
one extreme point in the decision space. Moreover, there 
may be more than one extreme point in the decision space 
mapping onto the sa,me extreme point in the objective space. 
Since the efficient frontier in the objective space is 
piecewise linear curve, it is enough to restrict ourselves 
to extreme efficient points. Thus, it is advantageous to 
look for efficient extreme points in the objective space. 



CHAPTER III 


GEFERATIOH OR EPPICIEHT POINTS 

3.1 Introduction ; 

In the multicriteria problems es already mentioned, 
one is interested in the generation of efficient points. 

It was also stated that efficient points in the objective 
space will be more useful. Because of the polyhedral 
property of the payoff set, the set of all efficient solu- 
tion is determined the moment the extreme points of the 
polyhedron are known. Since the primary motivation to 
this research is the algorithm proposed by Aneja and Nair[l], 
this will be discussed first, followed by the algorithm 
proposed by us. 

3.2 Ane.i'a and Nairas Algorithm ; 

3.2.1 Introduction ; 

Aneja and Nair proposed the following algorithm 
for generating the entire set of efficient points in the 
criteria space. Though the authors have concentrated on the 
determination of efficient extreme points for the bi-criteria 
transportation problem, th^ claim that their algorithm is 
equally applicable to bi-criteria LP problems. The algori- 
thm initially tests for the existence of a superior solution- 
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If such a solution does not exist, 'the linear progracuning 
pro'blem is solved using objective functions which are 
positively weighted average of the objective functions. 
This process is carried out till all the efficient extreme 
points are determined. 


^ 1-^1 


3.2,2 Algorithm ; 

Step 0 ; Find = Min (z^ j x^S) and z^^^ = Min 

and xG S). Record Z 2 ^^)and set k = 1. Similarly, 

find z^^^ = Min (z^^ x?rS) and z^^^ = Min (z^ | z^ = 
and S). If z^^^ ) , stop. Else 

record z^^b and set k = k+1. Define sets L =:|(1,2)^ 

and B = 0 , and go to Step 1. 

Step 1 ; Choose an element (r, s) S- 1 and set 

i 

i 

let X be an optiml solution to the linear program. 

Min. 


a'=-' = 

14^) 

- 4 ^) 

- 

Cl 4 ^ — 

(4^) 

-4=') 


^(r, s) ^ s) 


( 1 ) 


S.t Ax -4 b 

X ^ 0 

If there are alternate optimal choose an optimal 
solution X for which Cx is minimum. 

Let z^ = Cx and Z 2 = Dx 
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If (z , Zp) is equal either to zl^^) or 

(^1 » ^2^ ^ B = E U|(r, s)^ and go to Step 2. 

Otherwise record (zi^\ zl^^)such that = z^ and 

( k ) 

z^ = Z 2 » and set k = k+1. 


I = ZV^(r, k), (k, s)l 

Step 2 ; Set I = L - ^(r, s)| . 

goto Step 1. 


and goto Step 
If 1 = 0 stop. 


2 . 

Otherwise 


Aneja and Nair have proved the following results. 

Theorem 3.1 ; A point z^^^ = (z^^\ is an efficient 

extreme point in the objective space, if and only if 
(k) 

z is recorded by the algorithm. 

Theorem 5.2 ; The algorithm is finite and requires exactly 
(2k-5) iterations when it has ki'zZ) efficient extreme 
points, 

5.2.5 Numerical Example ; ' ‘ 

For better understanding of the Aneja and Nair’s 
algorithm, a numerical example considered by them is dis- 
cussed. Table 5*1 contains the data for a balanced trans- 
portation problem with two cost criteria. The two numbers 
in each cell represent the cost coefficients. The number 
in the north west corner indicates numbers 

in south-east comer denotes d^^’s. 



Table 3.1- Data for the problem 


Supply 

8 



The solution is indicated in Table 3.2. Tho 
values given by the side denotes the objective function 
values. 

Table 3.2; Optimal solution with c- - cost 

J- J 

coefficients 
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Similarly = 167 (= Z2 ( 15ie corres- 
ponding solution is indicated in Table 3.3. The objective 
function values are indicated by the side. 

Table 3.3; Optimal solution with d^^ costs 
coefficients . 



(2 


(2 


208 

167 


Iteration 

No. 1 




Step 1; 

(r, 

s) 0 L = 

: (1, 2) 



4"'’ 

= 98 

(= 1265 

- I67l ) 


„{r. 

= 61 

(= (208 

- 143 / ) 

Min. 

4"’ 

Ox + 

(r, s) 

cl 2 

Dx 


S. t X ^ s (2) 

The Table 3.4 indicates the modified cost 
coefficients ^ =1,... mn 

where m and n are the size of the transportation problem. 
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Table 3.4 j Cost for transportation problem 
in Iteration ITo. 1. 

Supply 
8 
19 
17 

Demand 11 3 14 16 



Table 3.5” Solution for problem corresponding 
to Table 3.4. 


v,5; (3 


) 


1 ) (16 


= 156 


Solution to the problem (1) is shown in Table 5.5. 


156, Z 2 = Dx = 200 


Since (z. , z_) is not either equal to (zp^, zP^) or 


\ -CJ-. 9 

.( 2 ) 


(zp\ (zp^), we have, ^ 2 ^* 

Now, k = 3+1 = 4 and L = pi, 2), (1,3), (3,2)j. 


L = L 


- 1(1, 2)( = p 


1,3), (3 


, 2 ){. 


Iteration No. 2 


(r,-s) = (1, 3) 


Summary of further calculations arc shown in the 


form of Table 3.6. 
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Table 3.6; Bicriteria Solution in the Objective 
Spac G. 



A computer program was developed to solve the 
problems by Aneja and lair's method. The solution shown 
in Table 3.6 required sivun linear program to "be 
solved. 

In the following section, ttie proposed algorithm * 
for generating the efficient extreme points in the 
objective space is described. 

3 . 3 An Adjacent Extreme Point B^cc jc’ur o; 

3.3.1 Background ; 

In Chapter II, it was pointed out that the deter- 
mination of efficient extreme points in objective space 
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rabher than the decision space offers two important advan- 
tages in case of nulticrit cria problems. Algorithn 
described in section 3.2 generates all efficient extreme 
points in tho objective space. Even though the convexity 
of the efficient frontier is recognized by Anoja and Eair [1], 
it is not fully exploited. It was also montioned that 2k -3 
linear programs need to be solved if k (k 2) efficient 
extreme points exists. 

In this section, an algorithm is developed which 
attempts to secure an adjacent efficient extreme point in 
each pivot operation. As explained later, the algorithm 
may not always be successful in generating the adjacent 
efficient extreme points. However, this does not d.ntruct 
the utility of tho algorithm since the next extreme 
efficient point will be secured in few more pivot stops. 

Before we describe the actual algorithm, certain prelimi- 
nary results wnich lay the foundation for the development 
of the algorithm are stated. 

3.3.2 Preliminary Results ? 

Let us define tho original problem a.s 

P 0 : Minimi zo ^CX, DX ^ 

s.t. AX < b 

X ?0 

where C, and D are n-dimensional vectors representing tho 
cost coefficients, of the two objectives. X is a 
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n-dimensional vector of variables. A is m x n tsatrix 
representing t jclinological coefficients and b is m x 1 
vector of requirements. 

Consider the following single objective, mathe- 
matical programs 


PI I Minimise CX 
s.t. AX-< b 
X > 0 


P2 ; Minimise Df 

s.t. AX <:b 
X ^ 0 


let the optimal values of PI and P2 be and 

X 2 

respectively. Problems P^. and P - are defined as follows 


Minimise IK 
s.t. AXZ b 
CX<^ z. 


X*^ 0 


Minimise CX 
s.t. AX 

DX 4 Z 2 
X 0 


let X* and X 2 oe the optimal solutions for the 

problems P^^f. and P^* respectively. Sadagopan [e ] has 
1 2 

suggested a procedure for reducing the PO problem to 

Pg ( 1*2 ^ problem paranctrization of 2 ^( 2 ^). The 

1 2 

procedure utilizes the following theorems. 


Theorem 3.3: 


In the optimal solution to the program P-g 
<5: [z*, ] the constraint DX 4 -Z 2 be a 

ono i.e, if x solve P- then Dx = z^. 


where 

2 

binding 



Proof ; 

« 

Assume the contrary i.e. 

'Dx ^ D£^ ( 3 ) 

2 * being the absolute ninimui:i of C over S 

Two cases will bo considorod. 

Ca se 1 ; let Cx^ - Ox (5) 

(3) and (5) tahen together contradict the fact that 
solves P * sincG x is a better solution. 

± 2i 

Case 2 ; lot CxJ 4 Ox ( 6 ) 

Consider the line segment [x, x£] which lies in the 
convex set S. Over this line segment the linear functions 
C and 1 are non-increasing. Because of strict inequality 
in ( 6 ), there exists a. feasible solution x* on the line 
segment [x, x^] such tiiat 

Dx < Dx’ <i ' 

2 (7) 

Cx >Cx' > Cx* 

Inequality (7) indicates tliat x' is a feasible solution 

to P- with Cx' <Cx which contradicts the optimality of x. 

^2 


Hence Dx 
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Theorem 3.4 : 

X (£ S is efficient if and only if x solves P— 

^2 

where 6 [z*, ]. 

The sufficioncy and necessity of this theorem is 
discussed hclow. 

Sufficiency ; Let x solve P— (henceforth whenever not 

^2 _ 

mentioned Z 2 ^ ['^2* ^^1 ^ to he efficient, then 

it needs to he shown that there does not exist any x € S 
such that, 

Case (i) Dx Ite and Ox ^ Cx 

( 8 ) 

Case (ii) Cx -if Cx and Dx Dx 

Case (i) ; Lot there he an x (x <<• S) such that Dx < Dx = Z 2 * 

Solution X is a feasible solution to P- hut not optimal since 

2 

the theorem 3.3 states that every optimal solution x° 
satisfies Dx° = •* 2 » Hence Cx Cx. Therefore, x is not 
efficient. 

Case (ii ) ; Consider some x ch S such that Cx ^ Cx, where 
X solves P- . Since x is optimal f'.l ^x cannot he foasihle. 

Zp 

(otherwise x would have solved ^5"^) • resilLts in 

Dx > Dx. Consequently x is not efficient. 

Thus X is efficient. 
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L-jt 5 be the Sct of e.ll efficient solutions. 

Suppose X (c E. Obviously Cx < Ox* us otlici'^viso x* will 

dominate x. let Ox = z^. k03v~>^ thut x does not solve 

P— but some other x' ^S solves P~ . By theorem 3.3? 

^1 _ ^1 „ 

Ox = z^ = Ox' . Sine 3 x’ is optimnl to P- , Px' <" lx,. 

^1 

Howover, beco,use x belongs to E this implies'- that 


Dx’ > Dx. Honce Djc' 

Y 


Dx. 


ConsGouGntly e: solves 


z 


1 


Thus it is possible to generate the entire s.-t of 
efficient solution.'^' bv p^ir.ame trie ally solving . Bj 
theorem 3.3, the generated solution will hav.? specific 
levels of attninmont (Zp), of th'.- second crit.';rion. 
Hadley [5] suggests that for the parametric vrriation 


of the right hand side of the constr^’int sot? duo.! °l-~'plex 
algorithm can b3 used. 


3.3.3 Algorithm (Adjacent Extroms Point Procedure) ; 

Sten 0; Solve PI and P2. Lot ^ and Zp be the optim.al 

solutions to PI and P? respectively. Solve P * and P *. 

^1 ^2 

let xj and xj represent the corresponding optimal 
solutions. 

St on 1 ; Lot C and D b^ the ’slativc cost vectors of tlio 
objective function of problen PI whose cost vectors are C and 
D. (In the first iteroticn, if alternate optimum to PI exists 
then choose among tbo rlternate optimum, that solution which 
minimizes the second objective with cost vector D). Deterraine 
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nn ad,]- cent extreme point solution in which clie entering 
v-i-iable in xj_ such that 


= Min 

j 



for c - 4 0 and d . > 0 


/ 


( 1 ) 


(In c-ioo of ti :■ break thou .arbitrarily). 

Set k = k+l, >->ecord = CX , z^"^^ = DX 


If rx =-■ Dz 2 , stop. OtherwiSGj, go to '‘-^tep 1. 

Thooroin 3.5 s 

Algorithm described a,bovo gonerates all efficient 
oxtr eno points. 

Proof ^ As mentioned in section 3 . 3»25 the basic idea of 
this algorithm is to generate all efficient extreme 
point.'- by solving E- over Z 2 ^ • Only one 

constraint (Px< Z 2 ) is added to PI. Hence it can be treated 
as right hand side parn.metrization problem for the added 
constraint. Dual simplex algorithm can be used to do the 
parametrization. 

Lot B be the- basis n'.trix associated with th:- problem 
PI. The matrix A with objective function in cononical form 
is represented as follows; 


■°B °II 



] 

z 1 

B ¥ ^ 


- % 


1 j 


where B (m x m) and H (m x n-m) are basic and not basic 
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mtric s o-r a, C is also divided into Cg (1 x n) and 
(l X n—m) reprcD :;nting' the basic r>nd non— basic cost 
vocto ns respec fcivoly, Vith the addition of the constraint 
Dx 4 z 2’ constraints can be written as follows. 

= b 


0 # s 


Vb Vf + 1*2 = 


( 11 ) 


Let A' and b' be denoted as in (12) and (13) respec— 


cively. B represents the basis of P- . fho’^ 

^2 


A 


F 0 




( 12 ) 


b 


L^2 J 


(13) 




B 0 

Dg 1^ 


(14) 


Taking tho inverse of B [6], we h-’ve 


B 


*-l 


B~^ 0 


-DgB“^ 1 


(15) 


Thus 


# A' is 


a? follows 


— 


r~ 




- 

B”^ 0 


B 

F 0 


I 

B^^F 0 

-DgH-^ 1 


Dg 



0 

%-DgB-llI 1 _ 


( 16 ) 
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ivhero Dg (ixEi) -^nd 1 ^ 7(1 ^ n-n) nre basic and non-basic cost 
vectors of D respectively. 


To deteraino the right hand side range of the added 
constraint, the system B* b’ ^ 0 is to bo solved. It 
essentially corresponds to the condition 


B 


whore 


DgB"-^ 

(k) 


0 

1 


n 

I 

1 


(k) 




0 


(17) 




represents the range over which B— continues 

^2 


to represent the optimal basis. Hence tho range will be 

{ ^ 2 , ^ 2 ] , where Z 2 is tho upper limit and Z 2 is ike lower 
Ck) 

limit of by solving the system (17). 


The variable to leave the basis will be the slack 
v'^riablo corresponding to the last row and the variable to 
'nter the basis will be decided by the minimum ratio rule 
of the dual simplex method. Since %-Dg F precisely 
corresponds to the r elativc cost a,s3ocir.tod with the 
objective function with cost .coefficients of D vector, the 
rule to decide the incoming variable in Step 1 of the 
proposed algorithm corresponds to tho minimum ratio rule. 


Theorem 3.6 ; 

Algorithm described in Section 5.3.3 terminates in 
finite number of iterations. 
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— <3u"l siciplex,* no basis is evor nopoaiod. In 
fact the previous basis becomes inf easible to successive 
ones. The number of extreme points in the objective space 

is finite. Hence the algoritlim terminates in finite number 
of iterations. 

It should b-^ noted that there may exist tie for the 

entering V'^riable at some iterations, i.e., more than 

ono v-^riable has the minimum latio of — c ./d . where 

3D 3 

and d^ are dual costs of non-basic variables. The tie can 
bo bromen arbitrarily. In such cases, next adjacent 
extreme efficient point ma.y not be obtained in one pivot 
operation. In case degeneracy is. encountered at some 
iterations, the lexicographic rule of simplex method can 
be adopted to prevent cycling. 

3.3.4 Numerical Examnlo ; 

The balanced transportation problem which was given 
as an illustrated example in Section 3.2.3, is discussed 
to bring out clearly, the features of the proposed algorithm. 
Data for the proble--^ is given in Table 3.1. 

Solution ; 

Step 0 ; Solving PI and P2, we h'^ve z| = 143, = 167. 

The corresponding solutions appear in Tabic 3.2 and Ihblo 3,5. 
Iteration Ho. 1 ; 

Step 1 ; Starting with the solution given in Table 3.2, 
the dual costs c. . and d . . of the non-basic colls are 

X J X J 
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det.?rrnlned , They ?,Ta shown in Table 3.7. 

Table 3.7; Dual COSTS for the transportation 
probloms at iteration no. 1. 


1 

V . 

3 

2 

V . 

3 


C5) 

'• 

l' 3'^) 

5 

-10 

r r 

3 

-5 

© 

3 

1 

.. _ ~5 


5 

10! 

1 6 



1 

2 

4 

4 






0 0 


0 


2 -8 


( 1 ) 

1 

( 1 ) 

2 


143 

265 


4 


13 


9 


In Table 3.7, the numbers in the north west come; 
indicates the dual cost c. - and the numbers in tho south 
cast corner represents tho dual cost d- u. and v. 
around the table represents u^ and v. with respect to 


2 - .2 


c. . costs. u3 and v. are u. and v. with respect to d . . 
13 13^3 

costs. 


Cell (2,3) enters the basis in the- next iteration | 

since the ratio of - c. ./d . is minimum for this non-basic | 

i J 1 .J I 

cell. I 

low Ir = 0+1, = 143. ^2^^ = 265. | 

Since is not equal to 167 (Dxp go to Step 1. | 

1 

I 

I 

[ 

# [ 
I 

I 
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Iteration Fo . ? : 


Step 1 ; 
revisod . 
function 


Now tliG dual costs of non— basic v'^riables ars 
liifey appear in Table 3.8xjit:-, the objective 
voluos by the side. 


Table 3.8; Dual costs of non-bnsic cellgat Iteration IIo.2. 



Purtlier calculations are shown in Table 3.9. 


Table 3.9s Solution to the transportation problem using 
the proposed algorithm. 


Iteration 

No. 

/ 

Leaving 

Variable 

Entering 

Variable 

r*" 

SI. No. of 
extreme 
efficient 
point 
(10 


1 

(2,4) 

(2,3) 

1 

(143, 263) 

2 

(1,1) 

(1,3) 

2 

(156, 200)' 

3 

(3,3) 

(3,2) : 

3 

(176, 173) 

4 

(1,2) 

(2,2) 

^ i 

(186, 171) 

5 

— 


- 5 1 

(208, 167) 
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Ths set of efficient extreme points were generated 
in 5 iterations. 0n3.y two line-r programs were solved 
and the rest of th.e calculations were only pivot operations. 
Whereas Aneja and Eair method reauiros solving of seven 
linear programs to arrive '^t the same set of oxtrenG 
efficient points. 





CHAPTER IV 


RESULTS are DISCUS SI ORS 

Computer codes were developed for the proposed as 
well as the Anoja and Rair’ s algorithms. The progrnjnnes were 
written in EORTRAR IV and implemented in DEC 1090 computer 
system. Since "both the algorithms require solving of 
linear programs, p, standard linear program code from 
Internationa]. Mathematical and Statistical library was utilized. 
The listings of programmes for Aneja and Rair’s algorithm 

t ^ ' 

and the proposed algorithm are given in Appendix A and B 
respectively. The execution tames for the randomly gene- 
rated problems of various sizes were recorded for both the . 

algorithms along with the number of extreme efficient points ; 
generated. It needs to bo pointed out that the lexicographic ^ 
rule required to prevent c.ycling due to degeneracy was not 
incorporated in the computer code. This results in repeti- 
tive generation of certain points. However, very few points j 
exhibited this characteristics. The tie for entering variable | 

I 

in tho proposed ".Igorlthm is resolvod using the stmtegy of 
steeper reduction. This aeons that of all tho non-bnsio ; 

variables, that is chosen tc enter the basis for which 
the product S Is aaximuffl. Here S represents the level at : 

J 

whicli X.: sn'fcGX'B "tho "bcisis# 

j 



31 


• 1 G-enorp.tion of Bandom Problems f 

In order to gonornte fop.sible '^nd bounded linenr 
progrniJSy balanced tr^^nsportation problens wore generp,tcd. 

The cost coefficionts generated were in the interval 1 to S. 

The demand a,nd supply wore in the interval 1 to 99 for smaller 
problems and 1 to 15 for 'higger problems. A separate Oonputor 
programme was developed for the generation of random problems. 
The programme listing is given in Appendix C. 


The randomly generated ba.lanced transportation problem 
of size m X n, is converted into a linear program with 
m + n - 1 constraints a,nd mn variables. Por example, a 
4x4 balanced transporta,tion problem results in a linear 
progrom with7constra3-nto and 16 variables, 

4 . ? Comnutation Timo s 


Computational experience for 24 randomly generated 
balanced transportation problems was gathered for both 
the algorithms. Problems of four different sizes viz., 

1- x4, 5x5, 6x6 and 7x7 were generated. Further, six 
problems for each problem size were consider^'d. The compu- 
tational times for solving these problems (in case of both 
the algorithms) are given in Table 4.1. For each problem 
size, Table i.2 gives the average execution times for the 
generation of one efficient point for both the ?.lgcrithms. 
The l«st column of this table gives a ratio risasure to 
compare the computational efficiency of two, .algo rir 




► 
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Table 4.1' Exacution ticos by the proposed Algorithin 
for problems examined. 


Problem 

Size 

No. of 
efficient 
extreme 
points 

Time taken 
for g€;ne- 
ration of 
extreme effi- 
cient points 
by Aneja and 
Nair's algo, 
(time in 
10'-3sec . ) 

ITo. of points 
genez’^-ted by 
the proposed 
algorithm. 

Time t'^ken 
for genera- 
tion of 
efficient 
points 

( t ime in 
10-3scc.) 

4x4 

1 

8 

3397 

8 

394 


2 

2 

839 

2 

344 


3 

2 

2779 

2 

381 


4 

5 

2317 

6 

384 


5 

4 

2060 

4 

379 


6 

3 

1078 

3 

283 

5x5 

1 

2 

1645 

2 

665 

2 

4 

4008 

4 

770 


3 

5 

5260 

5 

810 


A 

6 

7351 

6 

925 


5 

7 

8465 

7 

869 


6 

9 

9220 

9 

808 

6x6 

1 

4 

8418 

4 

1642 

2 

4 

8181 

5 

1781 


3 

8 

17112 

8 

1866 


4 

9 

20723 

9 

2129 


5 

6 

10 

22635 

13 

2286 


13 

30853 

14 

2388 

7x7 

1 

6 

18155 

6 

2321 ■ 

2 

7 

28897 

7 

2337 


3 

4 

5 

6 

9 

9 

11 

30857 

29105 

38964 

11 

9 

12 

2588 

2326 

2686 


12 

42863 

14 

2957 


J3 


for the genor'^tion of one efficient extrone point. Of the 
24 prohlens considered, for three prohloms, there was 
repetitive generation of the sane point at sone iterations. 
For one proble'^ of size 6x6 there were three repetitions 
while for two problems of 7 x 7 there were two repeta *ons 
each. As mentioned earlier, the lexicographic rule to avoid 
cycling due to degeneracy was not incorporated in the 
computer code and this r esul ted in the repetitive generation 
of theso points. Further whenever a tie existed for the 
entering variable the next adjacent extreme efficient 
point Was not generated in one pivot operation. In case of 
tic a point was obtained on the boundary lino joining 
the present extreme point and the next extrone adjac*.no 
efficient point. This phenomenon was observed in one 
problem of size 4 x 4, two problems size 6 x 6 and one 
problem of problem size 7x7. Only one boundary point was 
generated in each of these prob_i,cTTio . 

The average execution tine for one efficient point 
is O-Iloulatsd to divicllilfc- the total exGOUtlorattoe for all 
the six pvohlems of particular sizo by the number of effi- 
cient extremo points cs determined by Aneja and Hair's 
-algorithm. The second and third column of Table 4.2 show 
the average execution time for one efficient point by Aneja 
and Hair's algorithm and the proposed algorithm. 



T'lblo 4 . 2: Averr.go execution tines for generation of 

one efficient point by both the nlgorithns. 


4o. of 
problems 

Size 

Average time 
for gener-'ition 
of one efficient 
point by Anej.a 
•and Fair's -iigo- 
rithn. 

(A) 

Average time 
for generation 
of one efficion"'' 
point by propo- 
sed algorithm 

(B) 

[ 

Ratio ; 

(A)/(B) 

1 

6 

4x4 

445.35 

77.32 

5.75 

6 

5x5 

1089.36 

146 . 88 

7.41 

6 

6x6 

22-19.38 

251.91 

8.92 

6 

7x7 

2777.48 

281.38 

9.87 






Prom Table 4.2, it is cleor that the proposed algorithm 
is computationally superior to Aneja and Fair's algorithm. Por 
the 4 X 4 problem, the later algorithm required approximately 
6 times more execution time. As the problem size incroasos, 
the proposed algorithm shows even better computational 
superiority over the Ancja and Fair’s algorithm. Aneja and 
Fair’s algorithm requires solving 2k-3 linear programs, where 
’\s. {'^2) efficient extreme point exists, whereas the proposed 
algorithm requires only two linear programs and a little 
more than k pivot operations. This explains the reason for 
the computational superiority of the proposed algorithm over 
the Aneja and Fair’s algorithm. 
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It should he- noted tho.t the above stated cosiputor 

tines are based on solving transportation problems using LP 

code. Further computationa.1 tins reduction can be 

envisaged if a good transportation code is utilized, 

instead of the general linea- programraing code. This is 

due to the fact, the transportation algorithm will not 

require the calculation of the dual cost for all the 

vari'^'bles after each pivot operation. Only the u-^ s 

and ' s along .^rith tho corresponding row and column 

3 , 

invoxvod in tho pivot operation need to be computed , 

again. • 

i 

4 . 5 Conclusions ; | 

i 

t 

In this thesis, an algorithm which generates efficient | 

points in the objective space is developed to solve bicri- j 

! 

terion linear programs. [ 

The algorithm exploits the convexity of the efficient I 
frontier in the objective space. The original hi criterion 
problem is reduced to a problem of parametrization using i 

the results of Sadagopan [8]. This resulted in the 
proposed algorithm to require solving of two linear programs | 
and a little over k pivot operations where k is tho number ^ 
of efficient points. On the other hand, Anoja and Uair's I 

algorithm needs the solution of 2k-3 linear programs as the , 

authors do not fully exploit the convexity of the efficient : 


frontier. 



The coapnrison of the average execution tine to 
genero,te one efficient extrema point by both the rilgorithns 
clearly indicated the conputo-tioml superiority of tho 
proposed algorithn over Ane.jc and IKir's algorithm. For 
instance's the smallest problem i.e., 4x4 balanced trazis- 
portation problem required only one sixth of the a.verage 
execution time of Anujo, and fair's .algorithm. The comput'^- 
tional superioritj^’ is even more significant when larger 
problems a,re considered. 

For vQTj large problems number of efficient points 
may be large. In these cases an interactive procedure could 
be developed, which incorporates the proposed algorithm, so 
that the search region of tho efficient frontier is reduced. 

This algorithm as such C'^ji not be extended to problems 
with more than two criteria. However, there is a need to 
develop algorithms for generation of efficient points in 
the objective space for multicritoria problems. 
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APPENDIX A 


Program Listing for the ProgosefJ Ai22^il52 

mmmrnwmm m fuf m m m mmm mimm m***mmm mmm **«-»¥«««?► ei *t 


00100 

00200 

00300 

00400 

00500 

00600 

00700 

00800 

00900 

01000 

oiloo 

01300 

01400 

01500 

01600 

01700 

01800 

01900 

02000 

02100 

02200 

02300 

02400 

02500 

02600 

02700 

02800 

02900 

03000 

03100 

03200 

03300 

03400 

03500 

03600 

03700 

03800 

03900 

04000 

04100 

04200 

04300 

04400 

04500 

04600 

04700 

04800 

04900 

05000 

05100 

05200 

05300 

05400 

05500 

05600 

05700 

05800 

05900 

06000 

06100 

ttltl 

06400 

06500 


C*** mm ROUTINE ^ 

common /Bl/M,N.MA(50),MB(503aCOSTl(50,50).ICOST2(50,f 
COMMON /B4/AC50,50),B(50),IA,RWC4000) ,IW(2000) 

COMMON /B11/BASISC50) .IBASISC50) 

COMMON /B13/IZCNT,IZ(50,2) 

COMMON /B6/XC50,503 
COMMON /TO/INPUT, lOUT 
COMMON /lOl/IOUTl , INPUT! 

COMMON /B25/ALPHA(5) 

LOGICAL END 

INTEGER BASIS, DBASIS.TIMEl ^ ^ 

DATA lNPUT,IOUT,IOlITi ,INPUTl/20,5,22f23/ 

DATA IA/50/ 

ENDS, FALSE. 

NINE^O 

OPEN CUNITsINPUT,FILE='TRANSP.DAT'> 

OPEN CUNITxIOUT.FILEk'TRANSP.OUT') 


OPEN (UNITsIOUTI , FILEs ' IIMING . REPj^ j ACCESSa * APPEND 5 


C 

c 

C 

c 


401 


C 


to 


OPEN (UNIT=INPUTl,FILEa*TIMlNG.OLDh 


1000 


TYPE 401 
FORMATC 'Give the number 
ACCEPT ♦.NOPROB 
00 1200 1 1 JJssl, NOPROB 
CALL RTIME(TIMEl) 

CALL ZERO 

CALL DREADCEND) ^ 

IFCEND3 GO TO 1000 

CALL SETUP 

CALL FIRST 

CALL SECOND 

TYPE ♦tIZCNT 

CALL fiMECTIMEl ,IZCNT) 

CALL PRINT(NINE) 

NINE»9 
GO TO 1 
STOP 

CALL UERTST 
CALL ZXILP 
END 


of problems') 


FIRST 


COMMON^/B1/M^N?MAC50) »WB(50) ,ICOSTIC50,50) ,ICOST2(50»k 
cMon /B4/A?56,50) .B(50) ,IA,RWC4000) ,I»(2000) 


COMMON /B6/X(50,50) 

COMMON /B13/IZCNT,IZ(50j2) 


DIMENSION PSOLCSO) ,DSOLC50) .COST (50) 

INTEGER BOW# COL 
MlsM-tN-lJMasMiN 
DO 10 I a 1,M2 
PS0LCI)=0,0 
DSOLCIlaO.O 

g 0ST(I)*0,0 
OWaCl-.l)/N+l ^ ^ 

COLaI-( (ROW-l)*N) 

COSTCI)»-ICOST1(ROW,CQL) 

8all zxIlpca,ia.b,cost,m2.o,mi.s.psol,dsol.rm,i«.ier) 

DO 15 I a J,M 


DO 15 



A*» I 


06600 
06700 
06800 
06900 
07000 
07100 
07200 
07300 
07400 
07500 
07600 
07700 
07800 
07900 
08000 
08100 
08200 
08300 
08400 
08500 
08600 
08700 
08800 
08900 
09000 
09100 
09200 
09300 
09400 
09500 
09600 
09700 
191 


15 


urn 


20 


CONTINUE 

BCM1+1)=»S 

CAliL ZX3LP(A,IA,B,COST,M2,0,Ml + l ,S,PSOU,DSOl. , RW . I « » lER 
Ckhh MATFOM(PSOL) „ , 

CALL XCHANGCPS0L,C0ST,X,M2,SIIM) 

CALL PRINT(2) 

TENP1=0,0 
TEMP2s0,0 
DO 20 I = l,M 
DO 20 J = 1,N 

TEMPlsTEMPl+ICOSTlCI.a)*XCI#J) 

TEMP2sTEMP2+ICOST2(I.J)4X(I#J) 

CONTINUE 

:i 


C 

C 


IZ(i,l)=TEMPl+0, 

IZ(l,2)sTEMP2+0. 

IZCNTsl 

CALL PRINTC13) 

RETURN 

END 




SECOND 

SUBROUTINE SECOND 


lOOOO 

10200 

10300 

10400 

10500 

10600 

10700 

10800 

10900 

UOOO 

11200 

11400 

11500 

11600 

11700 

11800 

1900 

2000 

2100 

2200 

2300 

2400 

2500 

2600 

2700 

2800 

2900 

3000 

3100 

3200 

3300 

3400 

3500 


too 


4 


505 

5 


10 


20 

30 



50,5 


/B13/IZCNT,IZC50,2) 
/B6/X(50,50) 


COMMON 
COMMON 

INTEGER ENTERfLEAVE 
LOGICAL STOP,UBCOND 
STOP*, FALSE, 

Ml*MfN-UM2sM*N - ^ 

CALL FDLEAV CENTER, LEA VE.UBCOND) 
IFCUBCOND) CALL ERRORCl),^ r 

CALL PIVOTCMI , M2 , A, B, BASIS, ENTER, LEAVE) 
CALL COMPUf 
GO TO 5 
return 

END 


SUBROUTINE PIVOT CM. N, A, B, BASIS, KpCOL,KEYROW) 
DIMENSION AC50,50) ,BC50) .BASIS (50) 


INTEGER BASIS ^ 

DENOMR* AC KEYROW , KEYCOL) 
DO 3 a = l.N 
IFCACKEYROg, 

DO 4 I *1,M 


,J).EQ.1,0)2,3 


IF[IjEQ;KfeYR0W)G0^T0„4 


ji r V I # w ) , NE ,0 ,0 ) GO TO 3 
CONTINUE 
ICHANG«J 
GO TO 5 
CONTINUE 

loRiATC’ SOME THING WRONG IN 51290') 

A CKEYROwl 5 ) *A CKEYROW . J ) /DENOMR 

B CkIyROW 5 »B C KEYROW ) /DENOMR 

CONflNUE 

DO 30 I ® 1,M 

IF C I. EO, KEYROW) GO TO 30 

kzfi5illfe?McSg).A(KE«0W,.. 

i C I ) *B C I ) -AC I .KEYCOL) »B CKEYROW) 
CONTINUE ^ „ 

DO 35 I *1,M 


A* 2 



13700 

13800 

13900 

14000 

14100 

14200 

14300 

14400 

14500 

14600 

14700 

14800 

14900 

15000 

15100 

15200 

15300 

15400 

15500 

15600 

15700 

15800 

15900 

16000 

16200 

16300 

16400 

16500 

16600 

16700 

16800 

16900 

17000 

17100 

im 

17400 

17500 

17600 

17700 

17800 

17900 

18000 

18100 

18200 

18300 

18400 

18500 

18600 

18700 

18800 

18900 

19000 

19100 

19200 

19300 

19800 

19900 

20000 

20100 

20200 

20300 

20400 

20500 


C f .I q jkIIrOW ) A ( KEYROW , KEYCOL ) = 1 * 0 
IFCBASIS(I),EQ,ICHANG)BASIS(I)*KEYCOL 
CONTINUE 
RETURN 
END 


SUBROliTIME XCHANG(PSOL,COSTjjX,N,SUM1 

COMMON /B 1 /NOSP,NUSHPT,MA( 50 ) ,MBC 50 J ,IC 0 ST 1 ( 50 , 50 ) ,ICOS! 


dimension RSDL(N) ,C0ST(N) ,XCS0,50) 

INTEGER R0W,CUL 
SUMsO.O 
DO 30 I = l.N 
R0 Wx(T-U/N 0SHRT +l 
C0U=1-C(R0W-»1)*N0SHRT) 

XCR0W,C0l.)sPS0l.{I) 

SUMsSUM + PSOLCD^COSTd) 

CONTINUE 

RETURN 

END 

COHSoN^/Bf/i^N?MAC 50 ),MBC 50 ),ICOSTl( 50 , 50 ),ICOST 2 ( 50,50 
COMMON /B 4 /A( 50 , 50 ),B( 50 },TAjRWC 4000 ),IW( 2000 ) 

COMMON /BU/BA5iS(50),IBASISC50) 

COMMON /B13/IZCNT,IZ(50,2) 

COMMON /B6/X(50,50) 

DO 10 I = 1,50 
MACI)=0?MBCI)=0 
B(I)»0 
BASIS(I)®0 


IBASIS(I)sO 
IZCl,l)=0jIZ(I,2)s0 
DO 10 J = 1,50 


xci#a)»o.o 

10 CONTINUE 

Ms 0 |M« 0 ?IZCNT«O 
DO 15 I ® 1,4000 
R«(I)sO.O 

15 CONTINUE ^ ; 

DO 20 J = 1,2000 I 

20 IK(J )=0 i 

RETURN 

ST M lil 

SUBROUTINE TIME(TIME 1 . IZCNT) | 

COMMON /B 25 /ALPHA( 5 ) 

COMMON /B 1 /mIn?MA( 50 )!!mb( 50 ),ICOST 1 ( 50 , 50 ),ICOST 2 ( 50 . 50 | 
DIMENSION ALPHACI( 6 ) 

integer TIMEI ,TIME2,TIMDJF : 

CAUL RTIME(TIME 2 ) | 

c S^SBHS?S¥f?T55fUwH»0(i).i=i,6) 

c'*” 5 rITECI&u|i;i|o (ALPHA0CI).Iji.6).IZCNT,TIMDir : 

150 FOBMAT( 6 A 5 , *ISiX,I 7 ,lX,l 7 , 2 X) 

RETURN 

END 

c*** MATFOM 




COMMON /B 4 /AC 50 , 50 ),B(& 0 ), 1 A,K" 1 «uuu^ 
COMMON /B 11 /BA 5 ISC 50 ),DBASIS( 50 ) 
DIMENSION BASEC 50 , 50 ),PSnLC 50 ) ,DS 0 LC 50 ) 


integer basis, IPSOLCSQ) ,DBASIS 
M1®M+N"*l t MzsM^N 
DO 10 I = l£M2,, 
IPS0L(I)=PS0LCl)+.005 


jssd 

DO 20 I a lpM2 


A“ 3 



20600 

20700 

20800 

20900 

21000 

21100 

21200 

21300 

21400 

21500 

21600 

21700 

21800 

21900 

22000 

22100 

22200 

22300 

22400 

22500 

22700 

22800 

22900 


20 


25 


35 

40 

45 


C 

C 

C 

C 

c 


3000 C*** 

23100 C 
23200 
23300 
23400 
23500 
23600 

23700 4 

2|800 

23900 5 

24000 
24100 
24200 

24300 6 

24400 
24500 
24600 

24700 7 

24800 
24900 
25000 
25100 
25200 
25300 
25400 
25500 
25600 
25700 
25800 
25900 
26000 
26100 
26200 
26300 
26400 
26500 

iim 
26800 
26900 
27000 
27100 

275§E 


8 

Q 


11 


13 


IFClPSOt<CI),l.E.O.O)GO TO 20 

iJSiJ + 1 

BAS1SCJ)=I 

CONTINUE 

IFCJ.EO.MDGO TO 40 

DO 35 I = 1,M2 

DO 25 J a 1,M1 

IFCBASIS(J),EQ,I)GO TO 35 

CONTINUE 

BAS1S{M1)=I 

GO TO 40 

CONTINUE 

DO 45 J = l.Ml 

DBASISCJ)=BASIS(J) 

BASE(I/a)aACI,BASISCJ)) 

CALL INVERT(BASE,M1) 

CALL MATMLTCBASE,A,B,M1,M23 

RETURN 

END 


SUBROUTINE INVERT (A,N) 

INTEGER PIVR.PIVC^ vi-kai 

DIMENSION A(50,50).PIVR(50>»PIVC(50),PIVE(50),YCSO) 

DO 4 I=liN 

PIVRCDaO 

PIVC(I)=0 

Ksl 

PXVE(K)=0.0 
DO 9 I«1,N 

IF (l!^EgtPIVR(L)) GO TO 9 
CONTINUE 
DO 8 J=1,N 
DO 7 L»1»K 

IF CJ.EQ.PIVC(L)) GO TO 8 
CONTINUE, 

A8=ABS(A(Ija)) 

CsABS(PIVB(K)L „„ „ 

IF (AB.LT.C) GO TO 8 
11*1 
J JaJ 

PIVE(K)sA(I,J) 

CONTINUE 

CONTINUE 

DO 13 I=1»N 
BaAlIfJO),,, __ 4 , 

IF CI.EO.II) GO TO 13 

CONTINUE 

CONTINUE 

If’^cLlEjN) go to 5 
DO 16 0«Ln 
DO 15 I=1#N 
IRaPIVRCI) 

tcs»pivcci5 


A- 4 



27600 

27700 

27800 

27900 

28000 

28100 

28200 

28300 

28400 

28500 

28600 

28700 

28800 

28900 

29000 

29100 

2^200 

29300 

29400 

29500 

29600 

29700 

29800 

29900 

30000 

30100 

30200 

30300 

30400 

30500 

30600 

30700 

30800 

30900 

31000 

31100 

31300 

31400 

31500 

31600 

31700 

31800 

31900 

32000 

32100 

32200 

32300 

32400 

32500 

32600 

32700 

32800 

32900 

33000 

im 

33300 
“3400 
i3500 
3600 
13700 
33800 
33900 
34000 
34100 
34200 
34300 
34400 
34S00 


15 

16 


19 

20 


C*** 

c 


20 

25 

30 


35 

45 


55 


5 

15 


25 


YCtC)aA(IR,J) 
DO 161=1, N 
A(I,0)=YC1) 

DO 20 1=1, N 

00 19 J=1.N 

JRaPIVRCJ) 

JCsPIVCCJ) 

y(JR3=A(I,JC) 

DO 20 J=l,N 

A(1,J)=Y(J3 

RETURN 

END 


SUBROUTINE MATMLTCBASE, A, B, Ml. M23^^ 

dimension BASE(50,50) ,M50,50) ,6(50) , A1 C50 ,50) ,B1 C50) 

DO 30 I = I, Ml 

DO 25 J = 1,M2 

SUM = 0,0 

DO 20 K = l,Ml , 

SUM=SUM+BASE(I,K)1'A(K,J) 

A1CI,J)=SUM 

CONTINUE 

DO 45 1 = 1,M1 


SUMaO.O 
DO 35 J » 

SUM » SUM 
B1 (n=SUM 
DO 55 I = 
B(I)=BICI) 

DO 55 a = 1,M2 


1 ,M1 

+ BASE(I,J)*B(J) 


1,M1 




R 
END 


»"5g|/i?§§';?0),B(50,,IA,gK(40OO,,I-t2D00) 
COMMON /611/BASISC50),IBA5XSC50) 


COMMON '/Bi/M,N.MA(50) ,MBC50) , ICOSTI (50 ,50) » IC0ST2 C50 , 50 
COMMON /B6/X(50,5O) 


COMMON /B13/IZCNT,IZ(50,2) 
DIMENSION PSOLC50),COST(50) 
INTEGER BASIS, DBASIS 
M1=M+N»UM2=M*N 


DO 5 I = 1,M2 
C0ST(I)»0,6 
PSOL(I)»0,0 

gg-JIT. i,Bi. 

ctED^XCHANG(PSOu![cOST,X,M2,SUM) 

ra;§?issi.o.o 

DO 25 I » 1 »M 

SUM? i suMiirf I » J) ncofi I ( I , J ) 

SUM2*SUM2 tX C I , J ) »IC0ST2 ( I , J ) 

?i?llgiSf,l)=SUM1.0.5 

cl5i’'^PRM?2)fCAlt£’pRINT(l4) 

RETURN 

END 


•*** fndemt *** 


SUBROUTINE FNDENTCENTER,STdP) 

?OmK /Bl/M!N?Mil!50),MB(50),lCOSTU50,50),_ICOST2C50,5( 


A- 5 



34600 

34700 

34800 

34900 

35000 

35100 

35200 

35300 

35400 

35500 

35600 

35700 

35800 

35900 

36000 

36100 

36300 

36400 

36500 

36600 

36700 

36800 

36900 

37000 

37100 

37200 

37300 

37400 

37500 

37600 

37700 

37800 

37900 

38000 

38100 

38200 

38300 

38400 

38500 

38600 

38700 

38800 

38900 

39000 

39100 

39200 

39300 

39400 

39500 

39600 

39700 

39800 

39900 

40000 

4010Q 

40200 

40300 

40400 

40500 

40600 

40700 

40800 

40900 

nm 

41200 

41300 

41500 


15 


20 


C 

c 


22 

25 

28 


30 

35 

36 

40 

104 


40 

m 

103 


COMMON /B4/AC50,50),B(50),IA,RWC4000),1W(2000) 
COMMON /Bll/BASIS(50),IBASIS(50j 
LOGICAL STOP 

DIMENSION 0vJBARC5O),CJPARt5O) 

INTEGER ENTER , ROW , COL , BASIS , DBAS! S 
STOPS, FALSE, 

MIsM+n-UM2=M»N 

EMTERsO 

rMINs9999, 

DO 5 lsl,M2 

CJBAR(I)s0,0 

DJBARCDsO.O 

DO 40 a5l,H2 

DO 15 I = l,Ml 

IF(J|EggBASISCI))GO TO 40 

DUMMYlsO,0?DUMMy2=0.0 
DO 20 I sljMl 
R0Ws(BASISCI)-1)/N+1 
C0LsBASIS(I)-CCR0W«n4N) , 

DIIMMY1 = DUMMY1 + IC0ST1CR0W.C0L)*A{I,J) 
DUMMY2sDyMMY2+ICOST2CROW,CQL)*ACI#J) 

CONTINUE 



ROW, COD' 


DJBARCJ)=IC0ST2CR0W,C0L) 
IFCCJBAR( ‘ “ 


GO TO 22 


DUMMYl 
DUMMY2 

(jj.GE.O.AND.DJBARCJ) .LT.O.O) 

fF(FMIN2c-Cj8AR(J)/DJBAR(J)))40,28,25 
ENTERsJ , . 

DO 35 I s i-Ml 

IFCAa,ENTER),LE.0,0)GO TO 30 

ratioisbcivaJi.otSSJ „ 

irCRATIOl.GE.XMlNDGO TO 30 
XMIMlsRATIOl 

IF(A(l,J),LE.0,0)GO TO 35 
RATI01=B(I)/ACI,J) - „ 
lF(RATiai.GE,XMIN2)G0 TO 35 
XMIN2aRATI01 

TFCXMInI. LT, 9999,0. AND, XMIN2,LT, 9999, 0)GO TO 36 
GO TO 40 

IF(CJBAi(J)*XMIN2,GE,DJBAR(ENTER)*XMINl)MINsJ 

ENTERS MIN 
CONTINUE 

WRITE! lOUT# 104) (BASIS (I), 1*1 »M1) 

FORMAT! "BASISJ%2qa5j 1X5) ^ 

WRITE(I0UT,101 ) !COBAR(J) ,0*1 »M2) 

WRITE(I0UT,102) !DJBAR(J) ,0*1 »M2) 

DO 40 I * 1,MJ 



C 

c*** 

c 


IF(EMTER,NE,0)RETURN 
STOP*, TRUE, 

RETURN 

END 




CDLEAV 

SUBROUTINE FDLEAV (ENTER, LEA VE ,UBCOND) ,.„r«e«.n/eA k, 

COMMON /Jl/W|Sl«H50>f^2(50),ICOSTl(50|50)jICOST2(50,5( 



A- 6 



41600 

41700 

41800 

41900 

42000 

42100 

42200 

42300 

42400 

42500 

42600 

42700 

42800 

42900 

43000 

43100 

43200 

43300 

43400 

43500 

43600 

43700 

43800 

43900 

44000 

44100 

44200 

44300 

44400 

44500 

44600 

44700 

44800 

44900 

45000 

45100 

45400 

45500 

45600 

45700 

45800 

45900 

46000 

46100 

46200 

46300 

46400 

46500 

46600 

46700 

46800 

46900 

47000 

47100 

47200 

47300 

47400 

47500 

47600 

47700 

47800 

47900 

48000 

48100 

48200 

48300 

48400 

48500 


10 


c 

c 


LOGICAL UBCON0 
XMIN=9999,OjUBCOND=,TRUE, 
M1=M4'N-1,*H2=M*N 
DO 10 I = l.Ml 

IFCACI, ENTER) ,LE,0,03GO TO 10 
RATI0=8 ( I ) /A C I , ENTER ) 
IFCRATIO,GE,XMIN)GO TO 
XftllNsRATIO 
UBCOND=. FALSE, 

LEAVEsI 
CONTINUE 
RETURN 
END 


10 


10 

101 


15 

102 

20 


c 


SUBROUTINE ERROR Cl) 

COMMON /lO/INPUT.IOUT 
GO TO (10,15,10) I 

WRITECI0UT,101) . 

FORMATCISX, 'INFEASIBLE SOLUTION OR WRONG DATA*//) 
GO TO 20 

WRITE(I0UT,102) , , 

FORMATCISX, 'UNBOUNDED SOLUTION*//) 

Xa-1 

ROOTaSORTCX) 

END 


SUBROUTINE DREAD 


112 


155 


156 

101 


SUBROUTINE DREADCEND) 

LOGICAL END 

COMMON /Bl/MjNxHiC50),MBC50),ICOSTl(50,50),ICQST2C50,50 
COMMON /B25 /WhA(5) 

READUftPul!! !2,ENDalOO) (ALPHAd) .1 = 1 ,5) , M.INTO,N 
F0RrtAT(5A5,l4,A2.I4) 

READ(INPUT,101) (MACD , 1=1 ,M) 

READ(INPUT,10l) (MBCJ),Jb1,N) 

DO 155 J = l,M 

READ ( INPUT, l6l ) (ICnSTl(J,K) ,K=1,N) 

-- - 1,M. 


113 


DO 156 J = 

READ (INPUT, 101 ) (IC0ST2CJ,K) ,K=1 ,N) 
FORMAT(40I3) 

WRITECtOUT,113)(ALPHA(I).I=lr5),M,INTO,N 

F0RMAT(//////5A5,I4,A2,I4///) 

DO 15 11=1,2 




N) 


102 

103 


FORiM?Aij30x!dHE^lfiPul^COSf'l«ATRir'*;i5;//iOX.'DESTN 


(ICOSTia,J),Jal,N) 
(IC05T2(1 ,J),J=1,N) 


C5) 

.X, 'SOURCE 
= 2.M 


1'9X,20I5) 


104 

10 

15 


2)ARITE(lbUT,104)I,(ICOST2(|,J),J=l.N) 


106 

107 


108 

too 


f6rmAt(e 

DO 10 I 
IFCII.EC 

Vokmu 

CONTINUE 

SrITE(I0UT.106) (MACD, 1=1. M) 

FORMAT Cl OX, 'SUPPLY' 7X,20I5) 

FORMAT hox. 'DEMAND' ,7X,20I5) 

WRITE(I0UT,107)(MB(J) ,J=1,M) 

F^rIatCIHdIox, 'S olutions by Anela and Nalr method'//) 

RETURN 

ENDS. TRUE, 

RETURN 

END 


A"* 7 



48600 

48700 

48800 

48900 

49000 

49100 

49200 

49300 

49400 

49500 

49600 

49700 

49800 

49900 

50000 

50100 

S0200 

50300 

50400 

50500 

50600 

50700 

50800 

50900 


:ir 


c*** 

c 


10 


15 


r 

c 


...00 

51200 

51300 

51400 

51500 

51600 

51700 

51800 

51900 

52000 

52100 

52400 

52500 

52600 

52700 

52800 

52900 

53000 

53100 

53200 

53300 

53400 

53500 

53600 

53700 

53800 

53900 

54000 

54100 

54200 

54300 

54400 

54500 

S4600 

54700 

54800 

54900 

55000 

55100 

55200 

55300 

55400 

55500 


10 

110 


20 


25 

120 

125 


30 


35 

130 

135 


40 

140 


50 

150 


245 

250 

255 


256 


60 

160 


SUBROUTINE SETUP 

COMMON /B1/M,N,MA(S0) ,MBC50 ) , ICOSTI (50 ,50 ) , ICOST2 C50, 50 
COMMON /B4/A C50, 150 >, BC 50) .IA,RWC 4000), IW( 2000) ,LPCNT 
MlsM+N-l 


M2=M»N 
DO 10 T 5 1,M 
BCI)=MACI) 

DO 10 a S 1,N 
IJ= CI-l)=f=N J 
ACI,IJ)=l,0 
DO 15 I 3 M+1,M1 
8CI)=MB(I-M) 

DO 15 J = 1,M 
IJsCJ-l)fN t (I-M) 
ACl ,103 = 1 ,0 
RETURN 
END 


SUBROUTINE PRINTCCODE) 

Common /il/izi(§0)]iz2(50 ) ,izcnt,il( 50,2) ,TiiCNT 

COMMON /B1 /M^NiMA ( 50) ,MB( SO ) , ICOSTl (50,50), IC0ST2 (50 , 501 

E8SS8S1 ^li?dlfe5?:i?5?LlX,RW(4000,,IW(2000,,I,PC«T : 

COMMON /B14/MINCQC(2),MINCQD(2) 

INTEGER CODE 



WRITE(IOUT,120) CX(l,J),a=l,N) ; 

DO 25 I * 2,M ! 

iSRITE(I0UT,125) (X(I, J) , J=1 ,N) 

CONTINUE . ^ 

F0RMAT(5X,'The solution matrix follows V/5X,20(F6. 1 ,2X) 
FOPMAT(5X,20(r6.1.2X)) 

wllTE(I0UT,l30)(X(l,J),J=i,N) 

WRITECiqUTrisS) (X(I. J) , 0*1 ,N) 

F0RMAT(lx, *The shippino matrix follows*/5X,20(F7,l,2X). 
F0RMAT(SX,20(F7,1,2X)) 

RETURN 


FORMAT^lOxI^his is not an efficient point*//) 
RETURN 


WRlfE(I0UT,150) 
FORMAT(50X, *THE 


'//) 


PROBLEM IS TERMINATED 

FORMAT (20x1* No ^ of ^LPs Solved* * ,I6/20X, *^*' 
WRITECI0«T,?50) are* *,I5,/40X,29(*- 


FORMATOOX, 'Points follow* */30x , — — 

WRITE(IOUT,25§)(IZl(I),IZ2(I),I=l,IZCNT) 

F0RMATCC29X,I5,5X,I5/55 

8llTlciOUT,l60)IZCNT.lZiaZCNT),IZCNT,ip(IZCNT) 

FORMAT(20xI *The efficient points are ZC*,I3, ,1) 
15X,*Z(’,13,S2)a *,I5) 


*//30x.* 


s *.I5 


A** 8 


55600 

55700 

55800 

55900 

56000 

56100 

56200 

56300 

56400 

56500 

56600 

56700 

56800 

56900 

57000 

57100 

57200 

57300 

57400 

57500 

57600 

57700 

57800 

57900 

58000 

58100 

58200 

58300 

58400 


70 

170 


800 

805 


810 

815 


820 

10005 


1000 


900 

1105 


c*** 

c 


RETURN 

fORMATClSX. ’initial extreme points foXlowV20X, 'Z(l ,15* 
|! Z(lr2)sSl5,’ Z(2.i)sM5,' Z(2,2) = M5//) 

WRITE(I0UT,805) 

FORMATCSOX, 'Solutions by parametrically varying r,h,sid 
WRITE(I0UT,120)(XCl,J),a=l,N) 

DD 810 I = 2iM 

NR1TE( TOUT, 125 )(X(I,J),J=l,N) 

CONTINUE 

IF(IZCNT=s2)GO TO 820 


WRITECTOUT,815) IZl (IZCNT) .IZ2CIZCNT) , - , 

iri-i»Marnronif.7„jj^ C( I , J)*XCI , J) ’ ,F10 .4 , 'CORRESPNG DCI,J)*' 


FDRMATC20X, 

1,5X,F10,4) 

RETURN 

WRITBCIOUT, 10005) . 

FORMATC50X,’The primal is Infeasbible *//50x, 'The proble 
2ted and next problem pursued') 

RETURN 

WRITE(IOUT,1105) 

RETURN 

WRITEC|0UT,1105) 

F0RMATC50X,’ WHOLE RANGE COVERED' 

RETURN 
END 


) 


A- R 



APPENDIX B 

imm mtcn Aisssiisi 


00100 

00200 

00300 

00400 

00500 

00600 

00700 

00800 

00900 

01000 

01100 

01200 

01400 

01500 

01600 

91700 

01800 

01900 

02000 

oiJoS 

02300 

02400 

02500 

02600 

02700 

02800 

02900 

03000 

03100 

03200 

03300 

03400 

03500 

03600 

03700 

03800 

03900 

04000 

otIoS 

04300 

04400 

nin 

04700 
04800 
04900 
05000 
05100 
05200 
05: 

05^ 
05500 
05600 

81188 

05900 

06000 

81188 

06300 

06400 


C 

c 

c 

c 


c 

c 


PROGRAM MAIN ROUTINE 

**** VARABLE DECLARATION ^ 

COMMON /B 1 /M,N,MAC 50 ),HB( 50 ).ICOST 1 ( 50 , 50 ),ICOST 2 C 50 , 50 ) 
COMMON /B 4 /A? 50 , 150 J, 6 ( 50 ),IA,R»#( 4000 ),IW( 2000 ),LPCNT 
COMMON /R 5 /IZ 1 (SO),IZ 2 ( 50 ),I 2 CNT,IL( 50 . 2 ),ILCNT 
COMMON /IO/INPUT,IOUT 
COMMON /B 25 /ALPHA( 5 ) 

COMMON /lOl/IOUTl 

INTEGER TIMEl.TIME 2 ,TIMDir 

LOGICAL END, ILEMTY.INFCON.SUPRIRr STOP 

DATA 1 NPUT,IOUT,IOUT 1 / 20 , 21 , 22 / 

I A 3«50 

OPEN (UNITsINPUT.DEVICEsDSK, FILES *TRANSP. DAT') 
OPEN(UNITsIOUT,FILE='TRANSP,OUTp 

□PEN (UNITalOUT, FILES 'TRANSP. OUT SACCESSs'^PPEND*) 
OPENCUNlTsIOUTl , FILE® » TIMING, REP *, ACCESSs 'APPEND* ) 

1 CALL RTIME(TIMEl) 

CALL DREAD (END) 

IF(END).GO..T0..1500 


c ** 

c*** 

c 


CALL FEACHKCINFCON) ^ 

IFCINFCON) CALL ERRORCl) 

SET UP A MATRIX WITH EOUALITY CONSTRAINTS 

CALL SETUP 


C*** FIND INITIAL EXTREME POINTS 

^ CALL INITALCSUPRIR) 

IF(SUPRIR) 11/10^^ 

11 CALL TIMECTIMEI,!) 

CALL PRINT(l) 

GO TO 1 


C*** 

^\lo 


CHOOSE NEXT IL ELEMENT. IF IL IS EMPTY TAKE NEXT PROBLEM 

?ALL ?HOs|iI}LEMTY,IPOSTN) 

IF(ILEMTY) GO TO 1000 


^ CALL STEP21 (IP0STN,SUM1,SUM2,IR,IS,ST0P) 

IF(ST0P)G0 to 1 

lit * IF thI sotUTiaS'ii |o/eff||ci|ni do ho record 
C*** ELSE RECORD, DO THIS BY STEP22 

1000 

1500 ^JglE(UKITsINPUT) 

B- 1 



06500 
06600 
06700 
06800 
06900 
07000 
07100 
07200 
07300 
07400 
07500 
07600 
07700 
07800 
07900 
08000 
08100 
08200 
08300 
08400 
08500 
08600 
08700 
08800 
08900 
09000 
‘ ,00 
!00 
09300 
09400 
09500 

09800 
09900 
0000 
0100 
10200 
0300 
0400 
0500 
0600 
0700 
0800 
0900 
1000 

1300 
U400 

11700 

11800 

12100 

12200 

{?I88 

12500 

12600 

12800 

12900 

13000 


CLOSECUNIT»IOUT) 
CALL UERTST 
CALL ZXILP 
E^D 


SUBROUTINE ZERO 


COMMON /BI/M,N,MA(50),MB(SO),ICOS11(50,50),ICOST2C50,50) 
COMMON /B4/AC50,150) ,8(50) ,IA,RW(4000) ,iM(2aOO) ,LPCNT 
COMMON /B5/1Z1C50) ,1Z2(50),IZCNT,IL(50,2) ,IUCNT 
COMMON /B11/BAS1SC50),IBASISC50) 

COMMON /B13/IZCNTD.IZC50,2) 

COMMON /B6/X(50,505 
DO 10 I = 1,50 
MACI)=0?MB(i)s0 
BCI)»0 ^ „ 

BASIStDsO 

TBASISCI)=0 

IZ1CI)=0 

iz?iJi)»onzci, 2 )=o 
DO 10 J = 1,50 
Aa,J3s0,0 
X(I,J)*0,0 
10 CONTINUE 

Ma0?Na0;IZCNT«0 
IliCNTaO 
IZCNTDaO 
LPCNTsO , 

DO 15 I = 1,4000 
R»i(I)s0,0 
15 CONTINUE ^ 

DO 20 J a 1,2000 
20 IW(J)a0 

RETURN 
END 
C 

SUBROUTINE^TIMECTIMEI rIZCNT) 

COMMON /B25/AIiPHAC5) 

COMMON /BI/mIn?MA(50),MB(S0),ICOST1CSO,50),IC0ST2C50,50) 

5iSl¥E^IQ5n*iio rnifHM 1) , 1 = 1 , 5 ) , H.N.TimlF, IZCNT 

150 F0RHAT(5A5,2X, I2» *X',I2»2X,I7,1X,I3,20X) 

RETURN 

END 

C*** SUBROUTINE DREAD ♦♦♦ 

C 


SUBROUTINE DREAD (END) 
IfOGICAL END 



112 


gEMMoSI N jKI } i?> . MB C 50 ) , ICOSIl ( 50 , 50 ) , 

COMMON /825/AltPHAC5) 

RisDCI!Sm®fl2,ENO.lOO)(At,PHA(I),I=l,S), M.INIO.M 

READCINPUTflOi) fMAn),lal,M3 
READClNPUTflOl) CMB( J) , J=l »N) 

DO 155 J a l,M 


IC0ST2(50,50) 
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13500 

13600 

13700 

13800 

13900 

14000 

141O0 

14200 

14300 

14400 

14600 

14700 

14800 

14900 

15000 

15100 

15200 

15300 

15400 

15500 

156O0 

15700 

15800 

15900 

16000 

16200 

16300 

6400 

1188 

6700 

6800 

6900 

7000 

7100 

7200 

7300 

7400 

7500 

7600 

7700 

7800 

7900 

8000 

8100 

8200 

8300 

8400 

8500 

,8600 

8700 

8800 

,8900 

.9000 

9100 

.9200 

.9300 

-9400 

19500 

19600 

19700 

19800 

19900 

20000 

20100 

20200 

20300 

20400 


155 

156 
101 

113 


102 

103 


104 

10 

15 

106 

107 


108 

too 


c 

c 


10 


IS 


c 

c*** 


c 


READaNPUTa013CICnSTl(J,K) ,K=1,N5 
DO 156 J = l.M 

READ C INPUT, 1 01 )CICnST2(J,K),K=l,N) 
FURMAT(40I35 

WRITE(IOUT,1133CALPHACI),Is1,5),M,INTO,N 

F0RMATC//////5A5,I4,A2,I4///) 

DO 15 11=1,2 

WRITE(IOUT,102) TI,(I,I=1,N) 



1»9X,20I5) 


3;5X,20I55 

FORMAKSX, 'SOURCE 
DO 10 X ® 3 M 

IFai,EQ72)^RITEClOUT,1043I,ClCgS|2a,J3,0=l,N) 
IFai.EQ,l)WRITECIOUT,104)I,ClCOSTia,J),U=«l,N3 
FORMAT a0X,I4,9X,20I5/) 

CONTINUE 

COMTINUE 

WRITECI0UT,1063 CMA(I),I=1,M) 

FORMATClOX, 'SUPPLY' 2X,20|5) 

FORMATC fox, 'DEMAND ' , 7X, 2015 ) 
WRITE(I0UT,1073CMB(J),J=1,N) 

WRXTE{IOUT,108) ^ ^ 

FORMATClHl ,30X, 'Solutions by Aneja and Nair method'//) 

RETURN 

ENDs.TRUE, 

RETURN 

END 


0E5fN,*>' 


COMHo”/Bf/SSA(50),HB(50),ICOSTl(50,50),ICOSI2(50,50) 
COMMON /B4/AC50, 150) ,B(503 ,IA,RW(4000) ,1M(2000 J ,LPCNT 
M1«M+N-1 
M[2aM*N 

DO 10 X * IfM 
B(I)=MA(I) 

DO 10 J ®^t»N 
IJ= (I*l)4N + J 
A(I,IJ)3l,0 
DO 15 I = M+lfHl 
B(I)=MB(I-M),^ 

DO 15 J 5* 1,M 
IU*(J-1)*N 4 CI-M) 

A(i,ia5si,o 

RETURN 

END 


COMMON^/Bf/M!^N!MA(lo)’jMB?50),lCOSTlC5g,SO),ICOST2C50,50) 

HSmBoS /b|/a?5O|E)!|(|o),iLrw(4Ooo),iwc|oo0^ 

COMMON /B5/IZ1(50),IZ2(50) ,IZCNT,IL(50,2) ,ILCNT 

COMMON /B6/X(50,50) ncnrri^sm 

DIMENSION COStHBO) fPSOLC 150) ,DSOL(150) 


Ml*M+N-l 
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20500 

20600 

20700 

20800 

20900 

21000 

21200 

21300 

21400 

21500 

21600 

21700 

21800 


IS 


•000 
22100 
22200 
22300 
22400 
22500 
22600 
22700 
22800 
22900 
23000 
23100 
23200 
23300 
23400 
23500 
23600 
23700 
23800 
23900 
24000 
24100 
24200 
24300 
24400 
24500 
24600 
24700 
24800 
24900 
25000 
25100 
25200 
25300 
25400 
25500 
25600 
25700 
25800 
25900 
26000 
26100 
26200 
26300 
26400 
26500 
26600 
26700 
26800 
26900 
27000 

27300 

27400 


25 


28 


too 


C 


M2=H*N 

DO 100 11=1,2 
DO 15 I 5l,M 
DO is 0 =l,N 
IJ=CI-1)*N+J 

IF(II.E6.2)COSTCIJ)=-ICOST2tI,J) 
IFClI.KQ.UCOSTaJJs-ICOSTl (I,J) 

CONTINUE 

CALL ZX3LPCA,IA,8,COST,M2,0,M1,S,PSOL,DSOL,RW,TW,IER) 
LPCNT— LPCnT+1 

CALL XCriANG(PSOL,COST,X,M2,TOTCOS) 

TW0»2 

CALL PRINTCTWO) 

IFCII,EO,1}IZ1(II)=-TOTCOS 

IF(1I,EQ,2)IZ2CII)=-T0TC0S 

BCMUi)=iTOTCOS 

DO 25 I = 1,M 

DO 25 J = 1,N 

IJ=(I-l)*fi +J 

A(Ml+l,IJ)aIC0ST2tI,J) 

IFClI,EQ,l)ACMl+l,IJ)=ICOSTl(I,a) 

C0ST(IJ|=-IC0ST2CI,J) 

IFjJII^E|,2) C05TCIJ)=-IC0ST1(I,J) 

LPCNT=LPCNT+1 

CALL ZX3LP(A,IA,B,C0ST,M2,0,Ml+l,S,P50L,DS0L,RW,IW,IER) 
CALL XCHANG(PS0L,C0ST,X,M2,T0TC0S) 

CONTINUE ^ , 

IF(II.EO,l)IZ2(II)*-TOTCOS 

IFCII,EO,2)IZ1(II5*''TOTCOS 

TtflO 5? 2 

CALL PRINTCTWO) 

CONTINUE 

ILCNT=1 

ilci,i5=i 

IL(1#2)«2 

IZCNT*2 

fFCIZlCiKEo!izi(2),AND.IZ2(l).EO,IZ2C2)) SUPRIR».TRUE. 

RETURN 

END 


SUBROUTINE STEP21 ( IPOSTN, SUMl ,SUM2 , IR, IS, STOP) 

COMMON /B4/Af56,p0),IC50),lil,RW(4000),IWc5000),LPCNf 
COMMON /B6/X(50,50) 

COMMON /B25/ALPHA(5) 

DIMENSION COSTC150),PSOL(150) ,DSOLC150) 
logical STOP 
INTEGER THREE ,TWO 
STOP * , FALSE# 

THREE = 3 
TW0«2 

IR»IL(IPOSTN,n 



10 


DO 10 I »1|M 
DO 10 

foSTCiai=lAl*ICOSTl(I,J) + IA2*IC0ST2(I,J) 
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27500 

27600 

27700 

27800 

27900 

28000 

28100 

28200 

28300 

28400 

28500 

28600 

28700 

28800 

28900 

29000 

29100 

29200 

29300 

29400 

29500 

29600 

IVX 

29900 

30000 

30100 

30200 

30300 

30400 

30500 

30600 

30700 

30800 

30900 

31000 

31100 

31200 

31300 

31400 

31500 

31600 

31700 

31800 

31900 

32000 

liXo 

32300 

32400 

WXo 

32700 

32800 

33000 

33100 

33200 

33300 

33400 

33500 

33600 

nm 

33900 

34000 

nm 


15 

25 

IS 

30 


50 

55 

601 


c 


c 25 
20 


CQST(15=-C0ST(I) 

CALL PRINT(THREE) 

W1=M+N"1 

M2=M*N 

LPCHTsLPCNT+l 

CALL ZX3LPCA,IA,B,C08T,M2,0,M1 ,TOfCOS,PSOL,DSOL,PW, IW, lER) 
DO 35 I s 
DO 35 J = iIn 
IJ=(T- l)»NtJ 

A(M+N,ia)sIAl»ICOSTl(I,J) + IA2*ICOST2 ( I r J) 

COST (IJ3=»IC0ST1CI,J) 

B(M+N)»-TDTCOS 

LPCNTsLPCNT+l 

CALL ZX3LPCA,IA,B,COST,M2,0,Ml+l,S,PSOL,DSOL,RW,IW,IER) 

TF(IER,EQ,U3)G6^fO 55 

SaMl=0,0;SUM2sO,0 

CALL XCHANG(PSOL»COST,X,B2,TOTCOS) 

CALL PRINTCTWO) 

00 50 I =l,M 
DO 50 J = 1 ,N 

suMi=suMi+icosTia,a)»x(i,j) 

SUM2»SUM2tIC05T2(I,J)»X(I,J) 

RETURN 

ST0P»,TRUE, 

TYPE 601,CALPHACIII) jrlll-ltS) , ^ 

F0RMATC//5A5,' IS INFEASIBLE*// ) 

CALL PRINTCH) 

RETURN 

END 


SUBROUTINE 5ETP22 (SUMl ,SUM2 , IR,IS) 

COMMON /B5/IZl(50) .IZ2(50),IZCNT»IL(50,2) fILCNT 

COMMON /B6/X(50,505 
REAL SUM1,SUM2, 

INTEGER FOUR, SIX 
F0UR=4;SIXs6 
ISUMlaSUMl f ISUM2«SyM2 ^ ^ 

DIFR1*ABS(5UMI-IZICIR)) 

DlFR2sABS(SUM2-lZ2aR)) 

OIFSi*ABS(SUMl-IZiCIS)) 

DIFS2bABS(SUM 2*IZ2(IS>) ,,^TPc^ f -k « nwn nrirc-j it i 

IF( (DlFRl ,LTtl .AND,DIFR2,LT.1 ) ,OR . (DIP Sl , LT, 1 .AND,DIFS2,LT, 1 

13)25,20 

r*ftT r rR I MT c iF^ooR 5 
WRITES THAT THIS POINT IS NOT EFFICIENT POINT 
RETURN 






IZCNTsIZCNT+1 

IZ1CIZCHT3*ISUM1 

IZ2aZCNT3slSUM2 

ILCNT*ILCNT+1 

IL vILCNT , 2 jbIZCNT 
ILCNTbILCNT+I 

ReJ^RDS THE EFFICIENT POINT 
CALL PRINX(23 ^ 

CALL PRINT(SIX) 

RETURN 

END 



34500 

34600 

34700 

34800 

34900 

10 

35000 

35100 

20 

35200 

35300 

35400 

35500 


35600 

c*** 

35700 

c 

35800 

35900 

36000 

36100 

36200 

36300 

36400 

10 

36500 

101 

36600 

36700 

15 

36800 

102 

36900 

37000 

37100 


|7200 


37300 

c 

37400 

37500 

37600 

37700 

37800 

37900 

38000 

38100 

38200 

c 


.1300 
3^400 
3R500 
38600 
38700 
38800 
38900 
39000 
39100 
39200 
39300 
39400 
39500 
39600 
39700 
39800 
39900 
40000 
40100 
40200 
40300 
40400 
40500 
40600 
40700 
40800 
40900 
41000 
41100 
41200 
41300 
41400 


10 

110 


20 


25 

120 

125 

30 


35 

130 

135 

40 

140 

50 

ISO 

245 

250 


SUBROUTINE SETP3 CIR,I5) 

COMMON /BS/IZl|50 ), 122(50 ),IZCNT,IL (50, 23, lUCNT 

IF(nL(ia5!ES,lB).ANO,(II,CI,2),EQ,IS))GO TO 20 

CONTINUE 

RETURN 

IL(I,13=0 

IL(1,2)=0 

RETURN 

END 


SUBROUTINE ERROR(I) 

COMMON /IO/INPUT,IOUT 
GO TO (10,15,10) I 
WR1TE(I0UT,101) , 

FORMATdSX, nNFEASIBIE SOLUTION OR WRONG DATAV/) 
RETURN 

FoiMlia^xI^UN^OUNDED SOLUTION'//) 

RETURN 

END 


SUBROUTINE PRINT(CO0E) 

M ^ii^fi!?6E?^Z2(50),XZCNT,IL(50,2),ILCNT 

COMMON /Bl/MjN“HAC50),MBC5n),ICOSTl(50,50),lCOST2(50,50) 

M /M/A(iiaio|:i?5§Lx,R-c4ooo)a«(2oo 

COMMON /814/M1NC0C(2) ,M1NC0D(23 
INTEGER CODE 

GO^TOa0?20;30!40?50, 60, 70, 800, 900. 1000, 1000, 1000, 1000.820)C 

wRiTKlio6T,iio5;i£in),iz2{u,izi(5)az2(2) , 

F0RMAT(15X, ’Superior solution exists /20X, ZCl,l)- 15, 2(1, 

1,15,' Z(2,1)»^,I5.' 2(2,2)=', 15//) 

5rItIcI 0UT,120) (X(l,J),J=l,M) 

SRITl(JoUT,f25) (X(I,J),J=1,N) 

CONTINUE . 

FORMATCSX, 'The,SPlu_-- 
F0RMATC5X,20CF6,1 ,2X3) 

Si I TEC I0UT,130)(XC1,J),J=1,N) 

8RITl(I0UT,i35) (XCI, J) ,J=1 ,N) 

^ORMAt'cSx# 'T he shipping matrix £ollows'/5X,20(F7»l,2X)//) 
FORMATf5X;20(F7. 1,2X3) 

RETURN 

FORMlTClOxJ^fhis Is not an efficient point'//) 

RETURN ^ 

FORMAT cioxJ^HE PROBLEM IS TERMINATED '//) 

solved: ',16/20X, //> 

fflcient points are: ' , IS,/40X,29C *i»,*)// 
B- 6 


FORMATCSX, 'The solution matrix tollows'//5X,20(F6,l ,2X3/) 
FORMAT"*' 

RETURN 




41500 

41600 

41700 

41800 

41900 

42000 

42100 

42200 

42300 

42400 

42500 

42600 

42700 

42800 

42900 

43000 

43100 

43200 

43500 

43600 

43700 

43800 

43900 

44000 

44100 

44200 

44300 

44400 

44500 

44600 

44700 

44800 

44900 

45000 

.5200 

45300 

45400 

45500 

45600 

45700 

45800 

45900 

46000 

46100 

46200 

46300 

46400 

46500 

46600 

46700 

46800 

46900 

47000 

47100 

47200 

47300 

47400 

47500 

47600 

47700 

47800 

47900 

48000 

48100 

48200 

48300 

48400 


255 


256 

60 

160 


7( 
170 


800 

8U5 


810 

815 


820 

10005 


1000 

900 

1105 




10 

15 

C*** 

c 


'V/30X, 'One%8x, *T 


,15, 


WRITECI0UT,255) 

F0RMATC30X, 'Points follows V30x , — - 

1 #0 ^ / / ) 

WRITECI0UT,256)C1Z1C1).IZ2(1) ,I=1,IECNT) 

FaRMATCC29X,15r5X,T5/35 

R F T il R M 

WRITE(IOyTa60)lZCWT,IZlCIZCNT3 ,IZCN7,IZ2{IZCNT) . 
FORMATtaOXi’The efficient points are ZC^,i3,’,l)3 ' 

ISX^'ZC ' ,13, ' ,23= ',153 

WRITECigUT, 1703121(13, XZ2C13,IZIC23,IZ2(2) 

F0RMATC15X. ’initial extreme points followV20X, 'Z(l ,13 = ''I5, 

1' Z(1,23=*,I5,' Z(2,1) = M5,' ZC2,23 = ',I5//3 
RETURN 

WRI TE ( lOUT , 805 3 

FOHMA'nsOX, 'Solutions by parametrically varying r,h,side'//3 
WRlTE(I0Ut,1203 CXC1,J3,J=1,N3 
DO 810 I = 2,M^ 

WRITE(IOUT,1253CXCT,J3,J=l,N3 

CONTINUE 

TFCIZCNTa=2300 TO 820 ^ , 

WRITE(IOUT,8153 IZl ( TZCNT 3 , IZ2 C IZCNT3 t^, 

F0RMAT(20X, 'MIN C ( I , J3 ♦Xt I , J) ' , FlO ,4 , 'CURRESPNG DCI , J3*X(I , J3 ' 
1,5X,F10,43 

RETURN 

WRITEdOUT, 100053 . - ^ 

FORMATCSOX, 'The primal is inf easbible'//50x, 'The problem is apor 

2ted and next problem pursued' 3 

RETURN 

WRITEdOUT, 11053 
RETURN 

FORMlIf^OxI^^WHiLE RANGE COVERED') 

RETURN 

END 


SUBROUTINE FEACHKdNFCON) 

cHmmON^/BI/M?N,MAC503 ,MB(503 ,IC05T1C50,50) ,ICOST2(50,50) 

MSHRT=0?MSUPb0 

INFC0N=.TRUE. 

DO 10 I = 1,M 
MSUPsMSUP + MACI) 

MSHHTsO 

IFCcSaCI :UT,0).OR,(MB(J3.V|*0).OR.dCOSTld,J).LT.03.OR. 


HSHRT ^ 

IFCMSHRT.NE.MSUP) 
INFCONs, FALSE, 
RETURN 
END 


GO TO 15 


505 FORMA 


igSS8Si^5Si/5gf ?iS) ' IM50 , 2 , . ilcw 

LOGICAL 

ILEMTIa, FALSE. , 

WRITEd6llT,50S) |ILCNT 



153 



48500 

48600 

48700 

48800 

48900 

49000 

49100 

49200 

49300 

49400 

49500 

49600 

49700 

49800 

49900 

50000 

50100 

50200 

50300 

50400 

50500 

50600 

50700 

50800 

50900 

51000 

51100 

51200 

51300 


C 

c*** 

c 


DO 100 I = 1, 

TFdLdiD.GT 

CONTINUE 

IliEMT7 = ,TRUE, 

RETURN 

IPQSTNsI 

RETURN 

END 


IdUCNT 
GT,0) GO 


SUBROUTINE XCHANG(PSOIi,COST,X fN.SUM, 

COMMON /Bl/NOSP,NOSHRTeMA(50) fM8(50. 

DIMENSION PSOU(N),COST(N),X(5O,50) 

INTEGER ROW, COD 

SUMsO.O 

DO 30 1 = dN 

IPxpSOLd)+0,5 

PSOLd)=IP 

ROWsd-13/NOSHRT +1 

COU=I- C CROW-1 ) ♦NOSHRT ) 

SUM=SUM 4 PS0l.d)4C0STd) 
X(ROW,COL)=PSOLd) 

CONTINUE 

RETURN 

END 


,ICOST1C50,50),ICOST2(50 
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APPENDIX C 


EiSSH!! iSK SSE«IS£i2D Si ESDSSE EjSBiSSS 


OOlOO 

00200 

00300 

00400 

00500 

00600 

00700 

00800 

00900 

oigoo 

oiioo 

01200 

01300 

01400 

01500 

01600 

01700 

01800 

01900 

02000 

02100 

02200 

02300 

02400 

02500 

02600 

02700 

02800 

02900 

03000 

03100 

03200 

03300 

03400 

03500 

03600 

03700 

03800 

03900 

04000 

04100 



05000 

05100 

05200 

05300 

05400 

05500 

05600 

05700 

05800 

05900 

06000 

06100 

06200 

06300 

06400 

06500 


c 

c 

C 


4 were generated, 
by 5 were generated, 
by 3 were generated, 


C 

c 

c 

c 


10 


12 


15 


16 

200 


205 


1949 

18 




31 


35 


40 


ifilR problems of 4 by 

flilR twenty problems of 5 

with ISEEp 2948513# twenty problems of 3 
IMPLICIT INTEGER (A-Zj 
WENSION SUPPLy(5O),DEMAND(50)#COSTl (50,50), COST2(50, 50) 

NOOFPR=10 

IOUT*21 

flllR SivfR Ire SR® first set of problems of 4 by 4 is 2948513 
iflSR second set of problems of 5 by 5 Is 4548511 

ISEED GIVEN I OR THE FOURTH SET OF PROBLEMS OF 7 av 7 IS 40205096 
TSBED given FOR THE FIFTH SET 0F''k0Bi:i^^^ 4y^6 15^149554176 

^5*^41 76 

NOTIME=l 

DO 95 I a 1, NOTIME 
IFCI.EQ.USIZE = 6 

if(i,ne,i5size=size42 

Na24CSIZE**2)+2»SIZE 
DU 90 OJJal.NOOFPR 
DO 10 II a 1,N 
R(II)a0,0 
ISEEDalSEED+l 
CALL GGUB(ISEED,N#R) 

TOTALsO? RINDEXaO 
DO IS II a 1-SIZE 
RINDEXaRINDE^tl 
IMTERaRCRlNDEXl^lOO.O 

PAUSE » FINISHED SUPPLY' 

DO 20 Ilal.SIZE 
IFClIaaSIZlfelGO TO 19 
RIND|XaRINDEX+l 

GO TO 16 



TOT ALsTOT WTk'P 

So'to”o •'^'’•‘’<^*‘'1' "OOIFXCDEMAND, II, TOTAL) 

P^^SE * Finished demand* 

DO SO Ilai.siZE 
DO 50 JJal.SIZE 
DO 45 J a 1,2,1 
RINDEX=RIND£X+1 
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CONTINUE 

CONTINUE 

WRITE(igUTaiO) ajj, SIZE, SIZE 

FORMAT?' Problem no, ',14,' of size ',13,' x ',13) 
r0RMATC40I3) 

WRITE! lOUT, 101 HSUPPI.y til), 11 = 1, SIZE) 

WRITEtlOUT, 101) (DEMAND ( II), II=l, SIZE) 

DO 82 JN=1,2 
DO 82 JP=l,SIZE 

IFCJM==1)WRITECIOUT,101)CCOST1CJP, II), 11=1, SIZE) 
IFCJN=a2)NRlTE(IOUT,101)CCOST2tJP,II),II=l,SIZE) 
CONTINUE 

PAUSE ' One problem is over' 

CONTINUE 

CONTINUE 

STOP 

END 

SUBROUTINE MODIFYCOEMAND, II , TOTAL) 

IMPLICIT INTEGER CA-Z) 

DIMENSION DEMANDtSO) 

MAX=DEMAN0(1) 

DO 10 1=1,11 

IFtDEHANDCI) ,GE,MAX) GO TO 5 

GO TO 10 

MAX=OEMAN0(I) 

INDEX = I 

CONTINUE 

0EMAND(INDEX)=1 

TEMP a MAX-0EMAND( INDEX) 

TOTAL=TOTAL + TEMP 

RETURN 

END 


C* 2 



A63022 

A63022 

Date Slip 


This book is to be returned on the 
Jate last stamped. 



CD 6.72.9 





